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ABSTRACT 


The  physics  of  the  Gruneisen  parameter  for  solids  is  investigated  from 
the  view  of  conceptual  and  physical  contributions*  The  relatively  simple 
physics  of  a  van  der  Waal's  gas  is  applied  to  solid  densities  yielding  a 
Gruneisen  T  ^  2. 

The  separate  electron  and  phonon  contributions  tor  metals  are  considered, 
and  it  is  sho\^Ti  that  except  fur  very  high  and  very  low  temperatures  the  phonon 
ccmtribution  dominates.  Procedures  for  evaluating  P  from  shock-Hugoniot 
data  for  isotropic  solids  are  treated,  and  the  necessity  for  applying  ultrasonic 
data  to  anisotropic  crystals  is  considered. 

The  role  of  the  more  exotic  thermal  excitations  of  solid  state  physics  in 
determining  T  is  mentioned. 
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I.  INTRODUCTION 


Consider  a  system  having  total  internal  energy  £,  pressure  P,  volume  V, 
and  temperature  T,  The  Gruneisen  parameter,  P  ,  for  that  system  is 
defined  by 

r  .  . 

where  the  derivative  is  to  be  evaluated  at  constant  volume.  As  given  by  the 
macroscopic  (operational)  definition  of  £q.  (1),  P  is  a  measure  of  the 
change  in  pressure  produced  by  a  change  in  ^stem  total  energy  under  the 
condition  of  constant  volume. 

As  evidenced  from  the  above  definition  T  is  a  thermo  (energy) -mechanical 
(pressure)  quantity,  and  as  such  is  expected  to  be  important  to  thermomechanical 
problems,  A  few  examples  of  such  problems  are:  shock  wave  effects  (e.g, , 
e^qplosively  driven  shock  waves  in  solids),  the  thermal  expansion  of  solids, 
and  the  rapid  heating  of  materials  due  to  their  absorbing  intense  pulses  of 
nuclear  radiation.  In  the  last  example  P  becomes  of  extreme  importance 
since,  for  radiation  pulse  duration  small  compared  to  relevant  times  for  acoustic 
transport,  the  induced  thermal  pressure  is  directly  proportional  to  P 

In  this  report  we  shall  consider  the  physics,  sometimes  on  a  microscopic 
scale,  responsible  for  the  Grimeisen  parameter  having  its  observed  values. 

There  is  good  reason  for  doing  so.  Presently,  in  the  United  States,  technology 
has  outpaced  science  to  the  point  where  design  engineers  frequently  find  them¬ 
selves  using  materials  which  have  unknown  rei^nses  to  the  environments  in 
which  they  are  to  be  placed.  This  is  especially  true  in  the  area  of  military 
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missile  technology  where  very  often  new  exotic  materials  are  required  to  survive 
radiation  environments  in  which  their  response  is  unknown*  Knowledge  of  the 
Griineisen  constant  physics  would  allow  an  engineer  to  make  valid  estimates  for 
those  materials  for  which  P  has  not  yet  been  measured. 

Understanding  the  physics  of  P  is  also  quite  important  to  the  understanding 
of  related  physics  problems.  One  such  problem  area  of  interest  is  that  of 
electromeciianical  shock  effects  in  metals.  If  one  analyzes  shock  propagation 
in  metalL'  in  terms  of  a  two  fluid  model  of  electrons  and  phonons,  then  the 
Gruneisen  parameter  of  (;ach  component  becomes  important.  Some  of  the  analysis 
presented  in  this  repor!  will  be  directed  towards  that  problem. 

So  that  the  uniniliatid  reader  may  be  eased  gently  into  the  problem  we  will 
begin  by  treating  the  relatively  simple  problem  of  the  Gruneisen  parameter  of  gases. 
After  that  we  will  consider  some  generalities  of  P  .in  particular  as  they  apply 
to  solids.  Finally,  calculations  will  be  made  of  the  Gruneisen  parameter  for  some 
real  systems  of  interest. 
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n.  GASES 


A.  Perfect  Gases. 

Consider  an  Ideal  monatomic  gas  with  equation  of  state  (EOS)  and  total 
energy  E  given  by 

PV  =  NbT,  (2) 


From  Eq8.  (1)  and  (3)  we  find 

&i  the  above  N  is  the  number  of  molecules  (atoms  for  a  monatomic  gas),  k  is 

-16 

the  Boltzmann  constant  (1. 38  x  10  ergs/* Kelvin),  and  ^  is  the  absolute 
temperature, 

5 

Looking  now  at  a  gas  made  up  of  diatomic  molecules  with  E  =  2 
2 

r  -  ^  •  Consequently,  a  given  amount  of  energy  introduced  into  a  diatomic  gas 
produces  less  of  a  pressure  increase  than  in  the  corresponding  monatomic  gas. 
This  result  is  easy  to  imderstand;  for  the  diatomic  case  some  of  the  energy  goes 
into  the  Internal  energy  associated  with  molecular  vibration. 

While  the  above  example  is  somewhat  trivial  there  is  infonoatlon  of  interest 
contained  in  it  because  the  diatomic  gas  has  a  hint  of  solid  state  behavior  in  it  by 
virtue  of  its  Internal  energy.  Suppose  that  we  were  able  to  put  the  energy  Into 
the  diatomic  gas  ^stem  in  a  time  small  compared  to  that  necessary  for  transitions 
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5 

between  the  translational  states  of  a  diatomic  molecule.  Since  £  =  2  NkTis 
3 

composed  of  2  kinetic  energy  associated  with  translational  motion  of  the 
center  of  mass  of  each  molecule,  and  Nkf  internal  energy  associated  with  vibra¬ 
tional  states,  the  P  for  the  system  would  appear  as  shown  in  Fig.  1  below. 


r  + 


I 


Figure  1.  Grimeisen  parameter 
as  a  function  of  time  for  beam 
pulse  width  time  t.  The  time  for 
translation  between  vibrational 
states  is  t^  . 


O 


-  t 


Fig.  1  shows  the  P  that  would  be  measured  as  a  function  of  varying  deposi¬ 
tion  time  for  the  system  of  diatomic  molecules.  It  is  to  be  understood  that  we 

-14 

are  only  talking  about  a  thought  experiment  here  since  t  10  sec.  As  an 

^  „  2 

alternative  to  the  fast  energy  pulse  the  same  results  (i.e. ,  P  -  J  for  a 
diatomic  gas  system)  can  be  achieved  by  imagining  an  omnipotent  hand  reaching 
in  and  maintaining  the  pre-radiation  separation  distance,  and  thus  the  pre¬ 
radiation  state,  between  the  atoms  of  each  molecule. 
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B.  Imperfect  Gasee. 


A  gas  iB  considered  to  be  imperfect  when  interactions  between  the  gas  mole¬ 
cules  must  be  considered  In  order  to  explain  experimental  results.  One  of  the 
simplest  models  of  an  imperfect  gas  is  that  given  by  the  van  der  Waals'  EOS 


Pv>  NhT(.+  ^- 


(5) 


(6) 


where  a  and  b  are  constants  dependent  upon  the  gas  under  observation*.  The 

constant  b  is  generally  considered  as  being  proportional  to  the  volume  occupied 

by  a  molecule  (and  thus  unavailable  to  other  molecules),  while  a  is  proportional 

1 

to  the  interaction  potential  energy  between  the  molecules  . 

Combining  £qs.  (5)  and  (6),  while  neglecting  terms  of  higher  order  than  those 
Linear  in  a  and  b,  gives 


(7) 


Applying  Eq.  (1)  to  £q.  (7)  gives  for  the  Grvmeisen  constant  (i.  e. ,  energy 
derivative  at  constant  volume)  of  the  van  der  Waals'  gas 

£q.  (8)  contains  the  very  important  result  that  the  intermolecular  potential 

^Alternatively  a  and  b  may  be^consitered  as  making  up  the  second  virial  coeMcient, 
say  B2,  with 
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(other  than  the  approximately  hard  sphere  interaction  respansible  for  the 
volume  exclusion  constant  b)  does  not  appear  in  P  .  This  is  simply  understood; 
at  constant  volume  the  average  intermolecular  distance  remains  unchanged  by  the 
addition  of  energy,  and  thus  the  potential  energy  from  intermolecular  interactions 
remains  unchanged.  We  are  thus  led  to  expect  that  the  intermolecular  potential 
will  not  appear  in  the  constant  volume  Gruneisen  constant  of  a  solid.  Such  is 
indeed  the  case. 

Suprisingly  enough  £q.  (8)  can  be  used  to  estimate  the  T  of  a  liquid  or  a 
solid.  One  simple  model  for  b  has  b  =  4V^  ,  where  is  the  volume  of  a  molecule. 
Applying  Eq.  (8)  to  a  solid  or  liquid  with  half  of  its  available  volume  occupied  by 
the  molecules  (and  half  otherwise  empty),  with  the  simple  model  for  b  mentioned 
above,  yields 

which  is  a  suprisingly  acceptable  result;  with  the  exception  of  Zirconium,  metals 

2 

have  Gruneisen  constants  which  vary  between  one  and  three  at  room  temperature 
and  atmospheric  pressure. 

It  is  possible  to  understand  the  success  of  Eq.  (9),  as  applied  to  a  solid  or  a 

3 

liquid,  by  considering  the  shock  compression  Hugoniot  EOS  for  a  condensed 
medium.  For  volume  changes  small  enough  so  that  only  the  first  nonlinear 
term  need  be  considered  we  can  write 

or  i^on  rewriting 

frt  ^  C  3  (10b) 
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where  the  subscript  zero  denotes  the  undeformed  state,  the  subscript  H  denotes 
Hugoniot,  and  Ay,  are  characteristic  of  the  material  In  question.  ^  Is  the 
mass  density,  and  the  primed  quantities  are  the  constants  associated  with  the 
rearrangement  in  Eq.  <10b).  The  point  of  interest  is  that  Eqs.  (5)  and  (10b)  have 
the  same  functional  form  (A  is  a  function  of  temperature),  and  thus  the  proper 
interpretation  of  b  in  Eq.  (8)  should  lead  to  a  proper  result  —  which  it  does. 
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nL  GENERAL  THEORY 

There  are  two  ways  In  which  to  approach  the  problem  of  the  Gruneisen 
parameter*  One  method  involves  purely  macroscopic  thermomechanical 
considerations,  while  the  other  involves  a  detailed  statistical  mechanics  investi¬ 
gation  of  the  system  in  question.  We  begin  with  the  macroscopic  approach. 


A.  Thermomechanlcal  Cansiderations. 

We  start  with  the  definition  of  £q.  (1).  For  V  and  T  independent  variables  we 
have  at  constant  volume 


<1I) 


It  is  a  simple  problem  in  calculus  to  show  ,  by  alternatively  considering  V  as 
a  function  of  P  and  T,  that 


With  the  result 


(12) 


The  derivatives  of  £q*  (12)  each  relate  to  easily  performed  qua  si -static  (i.e* 
non  shock)  laboratory  experiments: 


Thermal  expansion  coefficient : 

'>r\ 


fis  Bulk  modulus  =  -  V  > 


vM 

(13a) 

(13b) 

e  a 

m. 

(13c) 

hi  terms  of  e(  ,  and  Eq.(12)  becomes 


r* 


(14) 
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la  the  absence  of  rate  dependent  effects  for  the  energy  deposition  problem 
(e,g, »  see  the  discussion  relating  to  Fig.  1)  the  quasi-static  result  of  £q.  (14) 
theoretically  should  give  the  same  P  as  that  of  the  dynamic  deposition  experi¬ 
ment.  In  practice  such  is  not  alwayti  trve  however^  as  the  energy  deposition 
experiment  measures  P  by  measuring  the  pressure  of  the  pvopagating  shock 
caused  by  that  deposition,  and  effects  such  as  attenuation  and  dispersion  can  easily 
mask  the  correct  induced  pressure  (and  V  )• 

There  is  another  macroscopic  approach  of  interest,  especially  since  it  is 

related  to  work  done  in  the  next  section.  It  happens  that,  to  a  good  approxlmatioa 
5,6 

,  the  thermal  properties  of  solid  can  be  represented  by  its  ''characteristic 

5  6 

termperature'',  @  .  In  that  approximation  the  internal  energy  can  be  shown  * 

to  be  of  the  form 

(15a) 

where  f  is  a  function  with  common  form  for  all  solids.  Using  the  first  law  of 
thermodynamics,  and  the  definition  of  the  free  energy,  F,  of  a  system 

F=E-TS^  (16) 

where  S  is  the  entropy  of  the  system,  one  can  easily  show 

+  (17) 

Eq.  (17)  says  that  F  and  E  have  the  same  functional  form 

(15b) 

where  again  g  is  a  function  with  common  for  r  for  all  solids. 
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Again  using  the  first  law  of  thermodynamics,  along  with  Eq.  (16),  leads 


e-mx- 

If  we  break  up  the  internal  energy  into  a  zero  temperature  contribution,  plus  a 
temperature  dependent  part,  then 

E(V^T>  »  (19a) 


£,C>^,T^-TS(nj>s5  F,  (Wb) 

Fc('4^  =  E,tV>j  (19c) 


where  Eq.  (19c)  follows  from  (19b)  and  the  third  law  of  thermodynamics 
(S  -»0  ae  T-»0). 


(20) 


where  the  second  part  of  Eq.  (20)  follows  from  the  characteristic  temperature 
approxlntatlon  as  given  by  Eq.  (15b). 

Operating  iq>on  Eq.  (15b)  it  is  a  simple  matter  to  show  that 


(21a) 


and  thus 


And  using 


i 

Hf)  ■ 


(21b) 


(22a) 
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in  Eq.  (17)  gives 


Substituting  Eq.  (22b)  into  Eq.  (20)  yields 


3V  \®  3V/  V 


(22b) 


(23) 


where  we  have  applied  Eq«  (22b)  to  the  temperature  dependent  cantribution  alone 
(subscript  1).  For  characteristic  temperature  models  such  as  those  of  Debye 
and  Einstein,  which  have  a  temperature  independent  0  ,  comparing  Eqs.  (23) 
and  (1)  gives 

n-  _  :iL  .  (24) 

e 

When  the  Debye  characteristic  0  is  used,  Eq.  (23)  becomes  the  Debye  EOS. 

By  using  the  fucctional  details  of  a  given  0  model,  and  ejqierimental 
Hugoniot  data  in  the  form  of  Eq.  (10a),  we  can  thus  use  Eq.  (24)  to  arrive  at 
numerical  values  for  P  .  Since  the  experimental  Hugoniot  data  need  not  be  from 
energy  deposition  experiments  (it  could  be  the  result  of  e]q}losive  Induced  shock 
measurements),  Eq.  (24)  turns  out  to  be  a  very  useful  resilt.  In  the  next  section 
we  shall  see  that 


P- 


(25) 


in  terms  of  the  coefficients  of  Eq.  (10a). 


B.  Statistical  Mechanics  Considerations. 

We  immediately  consider  a  system  made  iq>  of  two  weakly  interacting  subsystems. 
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As  such  the  formalism  Is  applicable  to  such  systems  as  fluid  filled  porous 

7 

solids  ,  mechauicai  mixtures  (e.g. »  KDNBF  with  diotomaceous  earth  -  an 

e3q>losive  mixture  sometimes  used  in  eiectro-eiqplosives  devices),  ard  for  at 

least  conceptual  purposes  to  electrons  and  phonons  in  solids. 

8 

Let  8  be  the  partition  function  for  the  system. 


2-  Z;  e 


rs 


2.  e 


Sfs  ) 


(26) 


where  E  is  a  total  energy  state  of  one  of  the  subsystems,  and  €  is  the  inter- 
r  rs 

action  energy  between  the  two  subsystems  when  one  subsystem  has  energy  E 

r 

and  the  other  has  energy  E  .  The  sum  is  overall  possible  energy  states  consistent 

s 

with  the  total  energy  of  the  system  being  constant.  The  interaction  is  assumed 

weak  in  that  E  and  E  are  taken  to  be  the  subsystem  energies  in  the  absence  of  any 
r  s 

Interaction.  The  mean  energy  and  entropy  of  the  system  are  given  in  terms  of  Z 
by 


E  -- 


'aJUi 


(27a) 


S-  k(JUi  +  ^E). 


(27b) 


Eqs.  (26)  and  (27)  are  sufficient  to  find  the  free  energy,  and  thus  the  pressure,  and 

finally  the  Gruneisen  parameter. 

-1 

Using  fi  =  (kX)  ,  and  fuiiher  defining  fi€  ^<1  for  a  weak  interaction,  Eq.  (26) 

rs 

becomes  to  first  order 


f.  (»•) 

_  — 

2  *  2a  P  ^  ^ts  ^  ^  J 

C  r  ^  fs  J  (28b) 
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where  the  subsystems  are  now  labeled  A  and  B,  and  is  the  partltioD  function 
of  subsystem  "A**  when  that  subi^stem  is  considered  to  be  an  isolated  system. 
Rewriting  Eq.  (2db) 


Thus 


and 


tni  = 

E*E*  — 


Performing  the  last  indicated  c^ration  in  Eq.  (30)  gives^^  ^ 


Similarly  allying  Eq.  (27b)  to  Eq.  (29b)  yields 


S-S.»S,+  ^  5  Sr. I 


£b  “ 


i 


—  -  ^1^3 


and  the  free  energy  takes  on  the  simple  form 


F=  F.+  F,+^ 


or 


F.  i-  F,  -  F, 


(29a) 


(32) 


AS  -> 


(33a) 


(33b) 


which  defines  .  Although  we  could  have  guessed  the  result  Eq.  (33b)  for  two 
weakly  interacting  sub^stems,  Eq.  (33a)  allows  one  to  calculate  Fy^g  when  one  has 
a  model  for  the  subsystems  and  their  interaction. 
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C.  General  Appllcatians. 

Ccmsider  the  case  (e,g. ,  a  mixture  of  weakly  interacting  gases)  when  the  sub¬ 
systems  can  be  considered  as  occupying  the  same  volume  •  Then 

P#  -  (34a) 

or 

fv  =  R.V  f  P,  V  -  V  ’  '*«*’> 

n  ■  +  r.i,  -V  ,  (34c) 

where  Eq.  (1)  has  been  used.  The  result  £q.  (34c)  is  quite  significant  in  that,  even 
for  ncm-intei'acting  subsystems,  it  says  that  the  Gnmeisen  parameter  is  a  weighted 
average  of  the  Gruneisen  parameters  of  the  sub^stems.  £q.  (34c)  can,  for 
example,  be  applied  to  problems  such  as  the  absorption  of  a  high  intensity  laser 

9 

beam  the  air-vapor  mixture  in  front  of  an  irradiated  solid  surface  . 

Let  us  now  apply  Eq.  (33b)  to  the  case  of  a  mechanical  mixtui'e,  the  subsystems 
(components)  occupy  volume  separately.  For  brevity  we  neglect  the  interaction 
energy  between  the  subsystems.  Thus 


F  =  F.  ^  (35a) 

V/  -  +  Vb  .  (35b) 


Eq.  (36a)  contains  the  realistic  assumption  that  does  not  depend  upon  Vq, 
similarly  for  F^. 


(36a) 

and 
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(36b) 


Pm  p  4.  p 

^  dV  ^  dv  ^ 


PV-  CE  -  PaV  ^  +  p,vi^ 


(360) 


and  then 


(36d) 


The  result  for  a  mixture  of  two  ncm-interacting  monatomic  gases  can  be  re¬ 
written,  from  £q.  (34c)  and  (2),  as 


an  expected  reeult.  £q.  (36d)  is  thus  reasonably  e3q)lained  as  Eq,  (37)  plus  a 

volume  effect.  We  note  that,  for  N  =  N*  N_,  V  ±  P  assumed,  we 

A  u  A  B 

have  P  varying  linearly  between  the  extremes  P  ,  P  . 

A  B 


For  the  q>ecial  case  of  components  A  aid  B  undergoing  equal  strain  (containing 
zero  strain  in  both  compcments  as  a  epecial  case),  the  solid  mixture  behaves  as 
the  gaseous  mixture  of  Eq.  (37).  For  equal  stx*ainB 


but 


(38a) 


(38b) 


■  4^  -  4^  «  . 

’  V  M, 


(38c) 
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Using  Eq.  (38c)  In  Eq.  (36d)  yields 

which  has  f  varying  linearly  betweoi  T  and  ^  since  E  =  E .  +  E  .  We  note 

A  B  A  B 

that  if  the  strains  are  equal  the  constant  volume  definition  of  Eq.  (1)  holds  even  for 
energy  deposition  pulse  duration  large  compared  to  the  time  for  acoustic  propagation 
across  a  component  particle;  if  the  strains  are  equal  they  will  be  zero  until  a  pres¬ 
sure  wave  due  to  inhomogeneous  deposition  reaches  the  local  observation  point. 

Thus  for  equal  strains,  which  contains  a  statement  concerning  elastic  constants, 

Eq.  (39)  appears  to  hold  for  energy  deposition  times  both  large  and  small  compared 
to  acoustic  propagation  times. 

There  is  another  q>ecial  case  which  can  be  treated  analytically.  If  the  component 
strains  are  not  equal  (thus  allowing  local  pressure  relief),  then  for  energy  deposi¬ 
tion  times  large  compared  to  the  time  for  acoustic  propagation  across  a  particle 
diameter 


f-  Pb 


and 


PV=  PftV  *  PaVa 


V, 


(40a) 

(40b) 

(40c) 


10 


where  indicates  the  "effective"  Grimeisen  constant  since  the  constant 

volume  conditions  of  Eq.  (1)  are  not  satistied.  Alternatively,  Eq.(40a)  can  be  used 
to  substitute  the  compressibility ,  K,  for  the  volume  terms  in  Eq.  (36d). 


-  V  JP 


(41) 
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When  Eqs.  (40a)  and  (41)  are  used  with  Eq.  (36d)  one  gets 

r*.  .  (42) 

'  'A  K£  ®  Ke 

ff  the  approximatians  of  Eqs.  (38c)  or  (40a)  do  not  hold,  then  Eq.  (36d)  must 

be  used  along  with  a  model  for  V  “  (V)  and  V_  =  V_(V). 

A  A  B  B 
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IV.  SOLIDS 


In  this  section  we  consider  the  individual  contributions  to  the  Gruneisen 
parameter  of  a  solid.  We  consider  only  metals  because 

a.  The  results  for  metals  contain  the  results  for  dielectrics  as  a  special 
case. 

b.  As  mentioned  in  the  introduction  we  have  a  strong  interest  in  electromech* 
anical  effects  in  metals. 

c.  For  the  most  part,  aerospace  materials  are  composed  of  metals, 
dieleotrics,  or  their  mixture  (e.g. ,  composites).  Further,  in  present  technology, 
aerospace  applications  stand  the  best  chance  of  seeing  an  environment  where  a 
knowledge  of  the  Grimeisen  parameter  is  important. 

We  break  up  the  problem  into  electron  and  phonon’*  contributions  while 
neglecting  the  volume  dependent  zero  degree  Kelvin  cohesive  energy  in  keeping 
with  the  constant  volume  definition  of  £q,  (1).  Thus 

FCVjT)-  *>  Sej>  (43) 

where  the  subscript  e  denotes  electrons,  p  denotes  phonons,  and  C^p(V,  T)  the 
contribution  to  the  temperature  dependent  part  of  the  free  energy  from  electron- 
phonon  interactions.  Even  though  the  interaction  term  appears  to  be  crucial  to  the 
electromechanical  problem,  for  the  Gruneisen  parameter  problem  one  gets  good 
agreement  with  experiment  by  considering  the  electrons  and  phonons  to  comprise 
non-interacting  gases  which  occupy  the  same  volume.  Thus  we  neglect  the  inter¬ 
action  term  and  note  that  Eq.  (34c)  applies. 

*A  phonon  is  a  quantized^init  of  lattice  vibration  in  a  normal  mode  system  where 
the  vibrations  behave  as  a  set  of  non-interacting  oscillators. 
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A.  Electrons. 


From  Zharkov  and  Kalinin  we  have  for  the  electron  system 


(44a) 

(44b) 

(44c) 

(44d) 

(44e) 


In  the  above  m  is  the  mass  of  an  electron,  N  is  the  number  of  electrons  in  the 

-27 

system,  and  h  is  Planck's  constant  (h  =  6.63  x  10  erg  seconds).  Comparing 

Eq.  (44e)  with  £q.  (3)  shows  that  the  electrons  behave  as  if  they  were  an  ideal  gas. 

11 

Eq.  (44e)  can  also  be  derived  on  the  more  general  grounds  of  thermodynamics 
with  a  Planck  distribution,  which  then  gives  Eq.  (3)  as  a  special  case. 

For  the  phonon  contribution  we  break  the  problem  up  into  low  and  high 
temperature  regimes.  The  reason  for  the  separation  is  the  radically  different 
functional  dependence  of  the  energy  and  entropy  In  the  two  regimes.  The  formulae 
for  both  regimes  is  from  reference  eleven. 

B,  PhoncMis-Low  Temperatures. 


t,W,T).V 
Sj.(.V,TjT  •  V 
/.  Fp(V,T)=- 


•0  ^ 

(45a) 

<4Sb) 

\s  (-fee) 

(45c) 

3t-  (-k  c) 
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vPpCx.T:..- V 


(45d) 

(45e) 


Li  arriving  at  Eq.  (45d)  we  have  treated  the  &ctor  of  V  in  Eq.  (45c)  as  a 
constant.  This  is  because  the  factor  of  V  tells  how  much  material  11^  present 
without  saying  anything  about  the  structure,  while  the  pressure  is  determined  by 
the  structure.  The  electron  subsystem  was  treated  in  the  same  way  in  that  the 
number  of  electrons,  N,  was  taken  to  be  constant.  hS.xit'Tt  ,  and  c  is  defined  by 


Reproduced  from 
uesi  available 


copy. 


(46) 


where  c^  is  the  longitudional  sound  velocity. 


and  c^  is  the  transverse  sound  velocity. 


C.  Phonons  -  High  Temperature. 

At  high  temperatures  the  sums  over  the  vibrational  states  become  greatly 
11 

simplified  with  the  result 


Ep  Cv^  T!)  =  3l4KTj 

SpCV,T:)T  *  -  Ut«  w  -h  1  j  ^ 


(47a) 

(47b) 


where  N  is  the  number  of  atoms  (for  brevity  we  consider  a  monatomic  lattice)  and 
9  is  some  mean  vibrational  frequency  associated  with  the  lattice  structure.  It  is 
obvious  from  Eqs.  (47)  that  only  the  entropy  contributes  to  the  pressure. 


(48a) 

(48b) 
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If 


we  now  use  the  relationship  from  lattice  dynamics 


4 


VA)  s  ^ 

where  ^  is  the  magnitude  of  the  state  vector  ^  ,  then  we  see  that  the  high  and 
low  temperature  regimes  yield  the  same  result 


(45e) 


A  few  words  are  in  order  regarding  the  difference  between  T  =  2/3  and  the 

e 

result  Eq.  (45e).  The  energy  of  an  individual  phonon  can  be  written 


6  =  ^  c  ^  (50) 

where  p  is  ^'momentum**  associated  with  the  phonon  state.  The  form  of  £q.  (50) 
is  radically  different  from  that  for  the  energy  of  a  gas  molecule  in  an  ideal  gas, 
namely, 


a.  tA 


(51) 


where  m  is  the  molecular  mass.  Thus  the  difference  between  T  and  P  is 

e  p 

explained  by  the  non-gas-like  behavior  of  the  phonon  system. 

From  the  discussion  of  section  n  of  this  report  we  can  assume  that  P  is  of 

P 

order  unity  (we  shall  actually  determine  it  later  in  this  section),  and  apply  Eq.  (34c) 

with  the  interaction  term  taken  as  zero  to  judge  the  relative  importance  of  the 

electron  and  phonon  contributions.  We  need  only  consider  the  ratio.  £  (V,  T)/£  (V,T). 

e  p 

The  ratio  is  calculated  below  for  T  =  300*K  (high  temperature  limit),  for  the  number 

of  electrons  equal  to  the  number  of  vibrating  atoms  (i.e. ,  a  metal),  and  an  electron 
21  3 

density  of  10  per  cm  .  The  low  temperature  i*atio  is  found  to  be  unity  at  10*K. 


23 


High  Temp. : 


EeC'^.T)  _ 

(52) 

EptViT) 

Ee  ^ 

(53) 

Ep(V,T) 

Low  Temp. : 


hi  Eqs.  (52)  and  (53)  we  have  taken  a  typical  solid  acoustic  velocity,  of  5  x  10  cm 
per  sec,  and  n  is  the  density  of  electrons  in  the  conduction  band. 

We  thus  see  that  a  t  room  temperature  only  the  phonon  contribution  is  of 
importance,  while  below  10*K  and  above  6,000*K  the  electron  and  phon«in  contri¬ 
bution  s  are  of  eqiial  magnitude. 

Let  us  now  calculate  from  experimental  shock  wave  data.  To  do  so  we  use 
£qs.  (10b),  (45e),  and 


-  fiET' 
"4 


(54) 


hi  using  £q.  (10b)  for  calculating  T  we  are  assuming  that  the  density  variation 

P 

of  is  primarily  due  to  phonon  effects.  Since  for  metals  the  electixm  ^*Bound'* 

5  7 

(i.  e. ,  Fermi)  velocity  is  of  order  10  cm  per  sec,  which  is  about  a  factor  of  20 

3 

larger  than  the  observed  propagation  velocity  for  a  weak  shock,  the  assumpticm 
would  appear  to  be  valid.  The  assumption  is  also  consistent  with  the  discussion 
accompanying  £q.  (52).  The  arithmetic  Is  straightforward: 


where 


5a.  (I 

Ak  \  ^  y 


(55a) 


(55b) 


(55c) 
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is  assumed,  and 


c  dv  “  c  av  t 


^  i?  Bm  ^  7 

’  5  1*'^  a;  Si‘J' 


(55d) 

(55e) 


'^r  -  T,  ’ 


(25) 


where  we  have  evaluated  Eq,  (55e)  at  V  =  V^. 

Thus  by  using  shock  wave  data  generated  by  non-energy  deposition  experiments 
(e.g. ,  flyer  experiments)!  we  arrive  at  a  Griinelsen  parameter  for  use  in  energy 
deposition  induced  shock  work. 

Table  I  below  compares  P  p  for  three  metals  from  Eq.  (25)  with  the  Gruneisen 
parameter  based  upon  the  thermodynamic  quantities  of  Eq.  (14). 


TABLE  I 

COMPARISON  OF  SHOCK  AND  THERMODYNAMIC  DATA 


Aluminum 

2.17^ 

2.09 

Copper 

2. 00 

1.98 

Lead 

2.78 

2. 46 

a« 

Hugonlot  data  from  Table  IV  of  reference  2. 

b. 

Values  from  L.  V.  ADtshuler  et  al,  Soviet  Phys.  J.  E.  T.  P.  573  (1960) 

The  agreement  shown  In  Table  I  is  rather  good,  and  is  t3rpical  of  the  agreement 
which  one  would  get  for  a  more  extensive  list  of  metals. 
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D,  Volume  Dependence, 

It  Is  also  possible  to  include  a  volume  dependence  for  P  in  the  calculation 
leading  to  £q.  (25),  One  need  simply  not  evaluate  Eq,  (55e)  at  V  =  and  for 
higher  volume  terms  simply  add  on  additional  higher  order  terms  to  £q.  (10b). 

The  volume  dependence  of  P  is  indeed  quite  importaiit  as  can  be  seen  from  the 
data  for  aluminum  (from  ref.  2)  given  below 

rtvi>  -»)+  ^5  (j  -  •  (56) 

The  volume  terms  of  £q,  (56)  contribute  twenty  percent  to  f*  (V)  for  f ”  1-2, 
and  the  result  would  even  be  more  impressive  if  the  second  and  third  terms  on  the 
right  of  £q.  (56)  were  not  of  opposite  sign. 

We  should  at  this  time  comment  on  the  physics  of  such  a  volume  dependence, 
namely  where  it  comes  from.  In  £q.  (8)  let  b  have  a  volume  dependence  (to  be 
explained  shortly),  and  expand  P  as  a  function  of  volume  for  a  small  volume 
change  about  an  initial  volume  V^.  Thus 

r(v.-»Av>.  P(0+  II  37l  J  K  ' 

There  are  two  contributions  to  £q.  (57).  The  first  involves  the  volume  change 

alone  (i.e. ,  AV/V  )  and  is  explained  by  the  observation  that  as  the  volume 
o 

increases  the  molecular  volume  becomes  less  important;  since  pressure  is  a 
measure  of  the  statistics  of  particle  interacticms  with  the  walls  of  an  imaginary 
container,  the  presence  of  other  particles  which  occupy  volume  effects  those 
statistics. 
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The  seccmd  contribution  to  £q«  (57)  is  proportional  to  (^^/dV),  and  says  that 
the  molecular  volume  is  dependent  upon  the  density  of  neighboring  molecules.  This 
is  understood  by  saying  that  where  the  volume  dependence  of  the  Gruneisen  para¬ 
meter  is  concerned »  an  exponential  repulsive  Interaction  is  more  appropriate  than 

a  hard  sphere  interaction.  Indeed,  by  using  such  an  exponential  repulsive  term 
12 

Sirdeshmukh  and  Rao  were  able  to  derive  P  (and  its  volume  dependence)  for 

some  cr^'^stals  of  the  fluorite  structure. 

Before  leaving  this  section  we  wish  to  mention  another  method  for  arriving  at 

the  volume  dependence  of  the  Gruneisen  parameter,  namely  the  Dugdale  -  MacDonald 
2  3 

formula  *  .  That  formula  is 


r>-^ 


I 


(58) 


2  3 

Eq.  (57)  is  derived  for  an  isotropic  solid  at  0*K  by  essentially  taking  the  third 
derivative  with  rei^ect  to  volume  of  the  lattice  potential  energy  (which  is  assumed 
to  be  only  a  function  of  the  atomic  coordinates).  Eq.  (58)  works  because,  if  is 
the  potential  energy,  and  if  the  potential  is  quadratic  in  the  lattice  parameter,  then 


1^4 

»  ■■ 


av 


(59) 


(ul/2ir )  being  the  vibrational  frequency  of  the  0*K  unoccupied  vibrational  states. 
The  success  of  Eq.  (58)  is  understood  as  the  right  hand  side  of  the  proportionality, 
Eq.  (59) ,  is  the  prime  ingredient  of  Eq.  (48b). 

3 

Equation  (58)  is  one  of  a  whole  class  of  relations  which  utilize  macroscopic 
PV  data  to  arrive  at  a  Grmeisen  constant.  All  such  relations,  Including  Eq.  (25) 


27 


and  the  procedure  of  this  sectioii  used  in  deriving  lip  assume  an  isotr(q)ic  solid* 

This  means  that  the  Hugoniot  data  used  in  arriving  at  P  must  be  in  a  range 

where  pressure  is  large  compared  to  the  yield  stress  (stress  for  the  onset  of 

plasticity).  This  is  because  most  of  the  Hugoniot  data  is  arrived  at  via  one- 

dimensional  strain  experiments  (e.g. ,  flyer  plate  with  planar  geometry),  with  the 

13 

result  that  solids  are  not  in  an  isotropic  state  in  the  region  of  the  yield  stress  . 

The  yield  stress  of  a  typical  metal  lies  between  a  few  kilobars  and  tens  of  kilobars. 

E,  Ultrasonic  Data. 

U  frequently  happens  that  Hugoniot  data  does  not  exist  for  a  material,  hi  such  a 
case  ultrasmiic  data  can  often  be  used  to  evaluate  P  .  What  is  needed  experimen¬ 
tally  is  the  usual  apparatus  for  measuring  elastic  ccmstants,  and  additional  equip¬ 
ment  which  allows  for  those  measurements  to  be  carried  out  as  a  function  of 
pressure  (a  maximum  pressure  of  a  few  kilobars  is  usually  sufficient).  The  theory 
is  relatively  simple.  If  ^  the  stress  tensor,  then 

where  ^  ,  Is  the  strain  tensor,  and  C  and  D  are  the  second  order  and  third  order 
kl  ’ 

5 

elastic  constants  respectively  .  For  a  cubic  crystal  undergoing  an  isotropic  volume 

change,  p  =  IF  =  tr  =  C" 

^  ’  11  22  33* 

Thus  the  constants  C  and  D  can  be  related  to  the  constants  A  and  B  of  Eq.  (10a). 

H  H 

The  D's  are  obtained  through  the  pressure  variation  of  the  sound  velocity.  Eqs. 

(60)  and  (61)  when  combined  are  of  the  same  form  as  Eq.  (lOb). 
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While  the  Hugcmlot  EOS  aj^roach  to  the  Gruneiseiii  canstant  is  valid  for 
xaotropic  materials  (e»g. ,  polycrystalliiie  metals)  it  is  not  valid  for  anisotropic 
single  crystals  This  is  because  the  Gruneisen  parameter  is  really  a  second 
rank  tensor  and  should  be  defined  by 


£q«  (1)  is  then  seen  to  be  Eq.  (62)  applied  to  isotropic  solids.  Since  the  C's  and 
D's  of  Eq.  (60)  give  all  of  the  necessary  information  for  evaluating  the  various  ele¬ 
ments  of  r  ,  the  ultrascaiic  method  is  the  siq;>erior  tool,  hideed  application  of  the 

hydrostatic  Hugoniot  method  to  an  anisotropic  crystal  such  as  Zirconium  (hexagonal 

15 

close  packed)  is  now  known  to  give  i  values  inconsistent  with  thermomechanical 
results  because  of  anisotropic  strain  effects. 

We  are  essentially  making  two  points  here.  First,  that  in  addition  to  the  usual 
Hugoniot  data,  shock  wave  physicists  have  relevant  ultrasonic  data  at  their  dii^sal. 
Second,  that  for  anisotropic  media,  only  the  ultrasonic  data  is  applicable. 
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V.  DISCUSSION  AND  SUMMARY 


The  main  Intent  of  this  report  has  been  to  give  the  shock  wave  physicist,  the 
prime  applied  user  of  the  Gruneisen  parameter,  an  understanding  of  the  physics 
involved  in  that  parameter  as  well  as  the  techniques  in  calculating  it  from  related 
data. 

In  this  closing  secticm  we  make  a  few  points  which  we  hope  will  serve  to  tie  up 
any  loose  strings.  We  also  mention  some  relatively  esoteric  physics  contributions 
to  which  at  this  time  although  interesting  in  their  own  right  are  relatively 
unimportant  in  applications. 

Within  the  main  text  of  this  report  we  failed  to  mention  the  connection  between 
the  characteristic  0  Grimeisen  parameter 


p. 

'  © 


> 


and  that  arrived  at  through  phonon  considerations 


V 

vu 


(24) 


(48b) 


The  ccnnecticn  is  that  the  characteristic  temperature  Q  allows  a  characteristic 
frequency,  W  /2ir  ,  to  be  defined  by 


(63) 


Thus  if  W  has  the  same  volume  dependence  as  m  the  two  approaches  will  give 
m 

the  same  result. 

At  high  temperatures  the  detailed  definition  of  W  for  a  monatomic  lattice  is 
11 

given 


w. 


(64) 
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where  o<  is  an  index  over  the  possible  vibrational  states  of  the  solid.  The 

3 

allowed  number  of  U]L<  increase  as  u«^  (a  relationship  known  as  the  density  of 

states  so  that  we  e:q)ect  w  «  as  is  associated  with  the  highest  possible 

energy  state  for  a  characteristic  0  model.  Ilius  we  are  led  to  e?q)ect  equality 

in  the  two  approaches  at  high  temperatures. 

At  low  temperatures  a  probability  factor  is  present  (relating  to  the  probabiltiy 

that  a  given  state  is  occupied  ii*  equivalent  of  Eq.  (63).  This  probability 

factor  serves  to  reduce  the  ratio  ui/lit^to  zero  as  the  temperature  goes  to  zero. 

Since  W  is  no  longer  close  to  Ul  they  do  not  necessarily  have  the  same  volume 

m 

dependence  -  and  indeed  for  many  materials  they  do  not  have  the  same  volume 
dependence. 

Germanium  offers  a  good  example  of  the  complications  introduced  by  low 

temperatures.  Experimentally  the  thermomechanical  Gruneisen  parameter,  Eq.  (14), 

14 

of  Germanium  is  negative  in  the  region  between  10*K  and  40* K.  Yet  Is 

determined  only  by  the  crystal  structure  and  remains  positive  in  that  same  region. 

That  low  temperatures  present  difficulties  should  not  be  suprising  as  we  have 

already  seen  that  electron  contributian  to  is  equal  to  the  phonon  contribution  at 

about  10*K.  Those  interested  in  further  details  of  the  electronic  contributions  to 

16 

are  referred  to  the  work  of  Wallace  . 

There  is  another  conceptual  approach  to  the  physics  of  the  Gruneisen  parameter 
which  is  very  much  worth  mentioning.  It  is  correct  to  think  of  phonons  and  electrons 
as  representing  thermally  excited  states  of  a  crystal  lattice  in  the  same  sense  as 
kinetic  energy  is  the  thermal  excitation  of  an  ideal  gas.  But  there  are  a  host  of 
other  possible  excited  states  in  a  crystal  lattice  which  are  more  or  less  important 
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depending  iqpon  temperature,  volume,  and  the  types  of  atoms  making  the  crystal 

in  question.  Some  of  these  thermal  excitations  are  listed  below. 

5 

MAGNONS  -  elementary  propagating  magnetic 

waves  in  paramagnetic  materials. 

5 

EXCnONS  -  cotq>led  electron-hole  pair  states. 

ROTONS^*^  -  elementary  rotational  states  in  a 

superfluid. 

ELECTRON-  5 

ELECTRON  PAIRS  -  the  elementary  excitations  responsible 
for  superconductivity. 

18 

HELICONS  -  an  electromagnetic  excitation  in  a 

solid  state  plasma. 

The  point  Is  that  each  type  of  thermal  excitation  represents  a  mechanism  for 
storing  energy.  Like  electrons  and  phonons  the  energy  dependence  of  the  partial 
pressures  contributed  by  those  more  exotic  excitations  must  be  considered  in 
evaluating  the  Gruneisen  parameter. 

The  effect  of  the  more  exotic  excitations  is  presently  not  a  serious  question  for 
the  applied  shock  wave  physicist.  However,  as  environments  change  and  previously 
strange  materials  are  introduced  into  technology  . 
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